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ABSTRACT 


The  equations  of  motion  of  a  rocket  are  derived  by 
applying  the  fundamental  definition  of  the  derivative  to 
Newton's  second  law  of  motion.  Three  independent  cases 
are  considered:  motion  of  the  missile  center  of  mass,  ro¬ 
tation  of  the  missile  about  a  transverse  axis  through  the 
center  of  mass,  and  rotation  of  the  missile  about  the  lon¬ 
gitudinal  axis.  The  second  cose  describes  the  motion  in 
pitch  (or  yaw),  and  the  third  case  describes  the  rotation 
(spin  or  roll)  of  the  missile  due  to  a  ring  of  small  jets 
placed  around  its  circumference,  The  equations  of  motion 
of  the  center  of  mass  are  then  modified  to  describe  the 
motion  of  a  satellite  moving  around  the  earth  in  a  nearly 
circular  orbit.  Finally,  a  method  is  developed  for  com¬ 
puting  the  approximate  impact  point  of  the  missile  by  al¬ 
gebraic  means. 


I.  INTRODUCTION 

The  purpose  of  this  Memorandum  is  to  derive  the  equations  of  motion  of  a  missile,  for  an 
oblatc-npheroidal  rotating  earth,  from  elementary  concepts  and  to  express  these  equations  in  a 
convenient  and  na  aral  coordinate  system.  Spherical  polar  coordinates  are  used  to  locate  the 
position  of  the  missile  center  of  mass  relative  to  the  earth.  It  is  assumed  that  the  earth  is  shaped 
like  a  perfect  sphere,  but  that  its  gravitational  field  Ib  that  corresponding  to  an  oblate  spheroid. 

The  motion  of  n  satellite  is  determined  directly  from  the  equations  of  motion  of  the  missile, 
the  former  being  considered  a  special  case  of  the  latter. 
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It.  MOTION  OF  CENTER  OF  MASS 


Consider  the  origin  in  Fig.  1  to  coincide  with  the  center  of  the  earth.  Let  the  set  of  axes 
(*»  z)H  be  fixed  in  inertial  soace,  and  let  the  set  £  ’(* y z  ')  be  rigidly  attached  to  the 
earth.  Consider  that  £  ’rotates  at  a  constant  angular  velocity  ft  with  respect  to  £. 

In  Fig.  1,  m j  represents  the  mass  of  a  quantity  of  exhaust  expelled  from  the  missile.  The 
center  of  mass  (C.  M.)  of  tftj  and  the  missile  is  located  at  Pj,  whereas  the  C.  M.  of  the  missile  is 
located  at  P2.  The  position  vectors’?  and ^5  are  drawn  from  Pj  to  mj  and  P2,  respectively.  It  is 
assumed  that  sxp.  0,  The  position  vaetorsTj,  R.  nndT are  drawn  from  the  origin  to  mj.  Pj, 
and  Pj,  respectively.  Let  t  »  tg  be  the  instant  of  time  just  before,  and  t  m  t0  +  A  the  instant  of 
time  just  after,  mj  is  expelled.  Here,  h  >  0  ia  an  infinitesimal  time.  Thus, 


pUp)  =  0 

(la) 

f*(tg)  ”  fl(tg) 

(lb) 

rftg  +  h)  =  R  (tg  +  h)  +  p  (tg  +  h) 

(lc) 

^1  (<o  +  A)  **  RUq  +  h)  +  8  (tg  +  h) 

(Id) 

«1  ”  m  U0)  -  m(tg  +  h ) 

(le) 

where  m(tg)  and  m(t0  +  h)  refer  to  the  mass  of  the  missile. 

Let  ^represent  the  total  linear  momentum  as  measured  in  £  .  Then 

■p(in)  =  m(tn) 


p(t0  +  /t)  =  m]  ~T|(t0+ A)+ m(l0+ A)  J-T Up+A) 


From  Eqs,  (1), 


rj  (f0  +  h)  --  r(tg  +  h)  t-"s(t0  +  h)  +  It) 


°Tho  nomenclature  used  throughout  this  Momornndunt  is  defined  in  Table  1. 
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Substituting  this  and  Eq.  (le)  into  the  expressions  for”p(t0)  and"fT(t0  +  A)  yields 


ru0  +  h)  -~pU0)  =  m(t0)  r(t0  +  h)  -  r  (t0)j 


-  m  (t0  +  A)J  iL  [T(<d  +  h)  -"p  (t0  +  A)J 


F  *  external  forces  acting  on  the  system 

p  (tf)  +  A)  ~?(<q) 


A.-rUQ  +  h)-.irnt0) 

-  «(<0)  lim  — - _lL_ 

A-*0 


m(t0  + A)-m(t0)  d  [W  .  ,.l 

- 2’T - -  •  [  *(<o  +  A)  -  p  (t0  +  A)J 


F  -  m(t0)  -iiT(t0)  -  m(t0)  i-  [T(<0)  -7(«0>] 


In  Eq.  (2),  d/dt  signifies  rate  of  change  in  £.  The  quantities  F,  r,  s,  and  p  can  be  ex¬ 
pressed  in  either  set  of  coordinates.  Since  it  is  most  likely  that  F  and  r  wilt  be  measured  in  £  *, 
it  will  prove  more  expedient  to  express  Eq.  (2)  in  terms  of  rate  of  change  in  £  i.e.,  in  terms 

of  d'/dt.  Now, 


dr  d'r  •?  — 

=  -r-  +  n  X  r 

dt  dt 
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and  similarly  for?and  p  .  Therefore,  Eq.  (2)  can  be  written 

F  -  m(ln)  [" -?L_.  r(t0)  +  2  fl  x  ~  r(t0 )  +  !Tx  iTx  r(t0)| 

'  0  U2  V  .*  J 

\ 

-  m(t0)  ["*"(<0)  "T (4o\j  +*^  x  s(<q)^  (3) 

since  0  »  a  constant  and  p(t0)  -  0,  Let 

—  [T(I0)  -*pl‘o>]  r.T 

where?  is  the  velocity  of  mj  with  respect  to  the  missile  (exhaust  velocity),  Now, 

1*21x71*0)1  '£  |7|  |7(i0)|.  For  s  typloal  missile,  j  s(f0)  |  £  20  ft,  Since  |1T|  =  7.27  x  10“® 
radians/'sec,  then  | S  x  $  (t0)  |  -  0.00145  ft/seo,  which  is  negligible  compared  with 
|?l 2000  ft/seo.  Henoe,  the  lest  term  in  Eq.  (3)  can  be  omitted.  Since  no  restrictions  are 
placed  on  the  time  t0  except  that  it  occur  before  mj  is  expelled,  and  since  is  being  continually 
expelled  (during  burning),  then 


7-  m(t)  [  i-  V  (t)  +  2  n  x  v(t)  +7x7x  r(f)]  -  M(t)  c  (4) 

where"T(t)  ■  (d‘/dt)~r{t).  Here  ,"7(0  and  F  are  measured  in  £  The  quantity  7 represents  the 
sum  of  all  external  forces  acting  on  the  missile  (excluding  vacuum*  ^rust),  and  the  quantity 
+  m(t)?equals  the  vacuum  thrust  of  the  rocket,  where  m  <  0.  The  quantities  -  2m  x? and 
-ml Tx  ft  x  r  are  the  Coriolis  and  centrifugal  forces,  respectively. 


Page  4 


Pog«  5 


Memoronrfum  Ho.  2Q-142 


Jtf  PnBultlon  Laboratory 


M 


li.  7(|.  ♦»-?<<!>)  , 

A-0  A 


—  p(tQ  +  A)  I" m(l0  +  A)  +  m(t0)l 
dt  L  J 


lim  ■<%+»— U»> 

A-.0  * 


4  (<0  +  A)  x  -^-4  (t0  +  A) 


+  Hm  V«7«(o  +  *) -/(*0)7V 

A— 0  J  A  . 

since  pU0)  m  0. 

£~<*o>  x  ^-p(*p)  [*•<*>]  -«(t0)T(t0)  x  j-  Td0)+  J-  /(<0)7(t0) 

if  ••  .. ' 

The  first  term  on  the  right  is  obviously  zero.  Now,  7(t0)  -7,  4(t0),  where?,  -  unit  vector.  Then 

^.T{*0)  -“a,  —  a(t0)  *  s(t0)  (y  xTj) 

whereT -  y<fix  .  This  is  illustrated  in  Fig.  2,  where  •>,  is  a  unit  vector  pointing  into  the  paper, 
and  Pj  is  a  unit  vector  lathe  direction  of  changing"?}.  Then^x"?,  -  Y$\  x*T,)  -  y<T,,  so  that 

* «0>  *  jj-  *<*0>  =  * «o> " l  x [i (t0)Ti  +  ys (t0)7j] 


=  ys2(t0)?j  x7, 

=  s2(‘o)  Y<t>\ 
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Therefore, 


(5) 


Assigning  a  value  to  s0  is  difficult,  since  s0  =  distance  from  P2  (at  time  t  -  tn)  to  mj,  and  since 
the  exact  location  of  mj  at  t  =  t0  is  not  known.  Let  s0  be  the  distance  from  P%  (at  t  ■>  t0)  to  the 
base  of  the  missile. 


IV.  ROTATION  ABOUT  THE  LONGITUDINAL  AXIS 

The  rapid  rotation  of  a  missile  about  its  longitudinal  axis  can  produce  the  same  stability 
in  flight  as  that  achieved  by  a  missile  employing  fin  deflection.  This  rotation  is  generated  by  two 
or  more  jets  situated  around  the  circumference  of  the  missile  at  the  roan  In  the  present  discussion, 
it  is  assumed  that,  if  f(  =  thrust  of  the  ith  jet,  then  51  f  ■  i  //  ■  0,  where  N  ■  total  number  of  jets. 

Figure  3  shows  the  rear  of  the  missile.  The  longitudinal  axis  passes  through  the  center 
of  the  circle  (representing  the  missile)  and  is  perpendicular  to  the  paper.  In  Fig.  3(a),  the  position 
vector  rt,  drawn  perpendicular  to  the  longitudinal  axis,  points  to  the  ith  jet  at  time  t  ■  (q.  Figure 
3(b)  shows  the  location  of  this  jet  at  time  t  «=  Iq  +  A.  The  quantity  is  the  moss  of  a  quantity  of 
exhaust  expelled  from  the  ith  jet,  and?^,  ff,  are  position  vectors  locating  m(,  as  shown.  The 
plane  determined  byT<(  s<(  is  perpendicular  to  the  axis.  From  Fig.  3, 


® i (*q)  "  ® 


•\ 


W  “  o) 


ft i (^  +  A)  m  r^(tQ  +  A)  +  $f(tp  +  A) 


>  (6) 
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Let  L  represent  the  total  angular  momentum.  Then 


L  (tg)  =  /  (Jq)  Cl)  (Jq) 


I«(t{)  +  A)  “  / Oq  +  A)  o>(<o  +  A)  +  \  R((t o  +  A)x  R((Iq 


+  A) 


*»1 


where  <u  «  angular  velocity  about  the  longitudinal  axis,  and  /  =  moment  of  inertia  about  this  axis. 
The  second  term  in  "?(t0  +  A)  is  the  angular  momentum  of  all  the  m(  about  the  origin.  If  it  is 
assumed  that  Wj  ■  mj  ■  =>  then 


it  .  .  n 

y  fy(to  +  A )xmi  R{(t0  +  A)  ■  mN  |?((t0  +  A)  x  r((t0  +  A) 
I«1  1-1 


+7,(I0  +  I.)  x  -^"T;(t0  +  A)  +"T^  (t0  +  A) 

at 


*  h)  +  A)*  +  A)] 


Now,  m(l0)  »  m(l0  +  A)  +  Nm/y,  where  m  refers  to  the  mass  of  the  missile.  But 


m  (l0  +  A)  -  m  (t0) 
lim  — . . . — - 

A-»0  * 


so  that 


m(t0)  =  m(t0)  »  -/V  lim 
at  .  _  A 


A-»0 


Paje  8 


4';?;.-,  ., 


X- 

~f::; 


MX: 

ffefrr- 

fife--. 

«iL. 

sjgwr 


if 

r>'  i*>' 


P- 

W&' 


JWaat,. 

.uJ 


vp 
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I.et  T  =  torques  taken  about  the  longitudinal  axis.  Then,  since  s^U0)  »  0, 


T  =  lim 
A»0 


L  (tQ  +  h)  —  L(t q) 


lim 

A-»0 


/(l0+  *)<?(*„  +  A)  —  /  (*0)  "£(t0) 


Nr  »  ,•* 

I,  mN  r-\  -  dr,  »  rfs, 

+  lim  _JL  \  r,  x  +  r,  x  —L 

A-0  h  L  *  „  *  1  * 


i  ■  1  u  / 

./ 


-  C  *7 
+  s^-rfT+s'x 


f  [«..)-(■„)]  -  ^  £>> 


i-i 


*[£T'(,°urH 


fiance, 


T-f['(,)"U)]  -iirZ'(^[r'“*r‘w] 


(7) 


But  -  r0^r,  where"r^ »  unit  vector  and  r  is  constant.  Also,  *«  Jj  <u(t),  where lies  along 
the  longitudinal  axis  (see  Fig.  4),  Hence,  (i/dt)  r(  "d)xrj=<i)r^  ,  where  &0  is  a  unit  vector 
in  the  direction  of  changing*^,  The  velocity  at  which  leaves  the  jet  iexhaust  velocity)  is 
equal  to  dT^/dt.  From  Fig.  4, 


ds. 

It 


dst 

It 


K) 


ve$), 
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Therefore, 

N 

2>«[ 

<-i 

+  -L T( (t)j  -  Nr(a>rZ0  -  veu0) 

so  that 

T- — 
it 

[fUaU]  «o  -  Ih(t)  r(ar~  vf)  3^ 

(8) 

Here,  <u(t),  rather  than  /<u(<)  it,  ia  of  ohiaf  interest.  The  torquea  T  are  highly  dependent  on  the 
particular  deaign  of  the  jeta. 


V.  EQUATIONS  OF  MOTION 

A.  Mlielle 

In  Fig.  5,  the  coordinate  axea  *,  y,  t  are  fixed  in  the  earth,  and  the  origin  coincides  with 
the  center  (of  mass)  of  the  earth.  The  quantities  i,  },  k  ore  unit  vectors  pointing  in  the  directions 
of  increasing  x,  y,  x,  respectively.  The  angular  velocity  vector  Q  of  the  earth  lies  alonjj;  the 
positive  y  axis,  as  shown.  The  position  coordinates  of  the  center  of  mass  (C.  M.)  of  the  missile 
are  r,  <£,  0,  defined  by  Eqs.  (9)  and  Fig.  5,  as  follows: 

*  =  r  cos  6  cos  <f>  v 

(y~ 

y  = r  </>  (  )  >  (9) 

V  -7T  S.  0  S  TT  J  I 
z  =  r  sin  6  cos  <j>  J 
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The  velocity  coordinates  of  the  C.  M.  are  v,  0,  a ,  defined  in  Fig.  6.  Here,?0,  <£0,  0O  are  unit 
vectors  pointing  in  the  directions  of  increasing  r,  <j>,  0,  respectively.  The  angle  between  the 
velocity  vector  7andT0  is  (77/2)  -S)  a  is  the  angle  between  <f>0  and  the  projection  of  7in  the 

plane  normal  to  rg.  From  Fig.  6, 


Now, 


but 


Hence, 


r  =  rn  r 


u  a  t»  sin  0  +  <f>Q  v  cos  0  cos  o  +  6q  v  cos  0  sin  o 


that 


d~  ~  ,  dro 
v  <=  —  «r0r+r 
dt  0  i it 


70  =*Tcos  6  cos  <js  +  T sin  ^4  +  T  sin  6  cos  tf> 

(jT0  »  i  ( -  cos  6  sin  <j>)  +  /  cos  <4  +  k  ( -  sin  6  sin  ^4) 
Tn  •»  f  (--ain  6)  +  k  cos  6 


-r-  -  <^o  I4  +  ^0  C0B  ^  & 

dt 


T  =  Tq  r  +  4>or  4>  +  "^o  r  008  ^  ^ 


(10a) 

(10b) 


r  (11) 


(12) 
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Equating  Eqa.  (10b)  and  (12)  yields 


r  «*  v  sin  0 


(13a) 


</>  =  —  coti  8  coa  a 
r 


(13b) 


A  v  «  8ln<7 

0  a  -  COS  O  - - 

r  cos  9 


(13c) 


From  Eq.  (10b)  and  Eqa.  (11), 


—  =<?0  (ti  coa  8  6+  v  ain  8  -  —  cos2  6)  +  0q  (- 1»  coa  8  ain  a  a  -  v  coa  a  sin  8  8 
at  r 


,  v2  v2  ,  , 

+  v  coa  8  coa  a  +  —  ain  0  coa  0  coa  a  +  —  coa2  ®  ain2  o  tan  <£) 


+  <?o  (v  coa  0  cob  a  a  -  v  ain  o  ain  0  6+  v  coa  8  sin  o  +  JL  ain  0  cos  8  sin  i 


„2 

—  cob2  8  coa  a  ain  a  tan  <£) 


(14) 


The  angular  velocity  of  the  earth  is  given  by 


Hence, 


0  =  ;'fl  =  Tq  ft  sin  <j>  +  ,^o  ft  cos  <j> 


ft  x  v  =  f0  ft  v  coa  0  cos  0  sin  a  +  <£0  ( -  ft  t>  coa  0  sin  4>  sin  o ) 


+  00  ft  v  (coa  0  ain  <j>  coa  a  -  sin  0  cos  <j>) 


(15) 
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n  x  0  x  r  =  r0  (-02r  cos2  4>)  +  4>o  sin  0  008  ^ 


(16) 


The  external  forces  F  acting  on  the  C.  M.  will,  in  general,  consist  of  the  following: 

1.  Aerodynamic  forces 

B  »  base  drag  effect 
D  »  drag 

•HV> 

L  «  lift  (excluding  rudder) 

/?*  -  lift  of  rudder 

2.  Nonaerodynamlc  forces 

•m 

1  ■  jet  vane  and/or  exhaust  obstruction  effeot 
mg  -  gravity 

It  Is  assumed  that  the  vacuum  thrust  /q  and  all  the  forces  listed  above,  excluding  gravity, 
lie  In  the  plane  determined  by  "u  andT0,  as  shown  In  Fig.  7.  Here,  x  S  »  /q  x  /  -  «x  D  »  Z,  ■  D  ■»  0, 
and  /V  is  normal  to  the  longitudinal  axis  of  the  missile.  From  Figs.  6  and  7, 


70  *  r0/0  sin  (6  +  a  +  ij)  +  cos  (®  +  a  +  »j)  cos  o  +  0O/O  cos  (0  +  a  +  ij)  sin  o 


-ft 


Jo~° 


y  U7) 


V 


L  cos  0  -  <^Q  L  sin  ®  cos  a  -  60  L  sin  0  sin  o 


a  =  -r0/y  cob  (©  +  a)  +  <^q  A  sin  (0  +  a)  cos  a  +  6q  N  sin  (0  +  a.)  sin  a 
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The  potential  function  for  an  oblate-apheroid  earth8  is 

.  GMm  F  R  R2  (  1  2  A  ^  lac  .  4  j  an  •  2  •  q\  I 

r,(4)  a  -  —  +  K,  —  ( - ainV)+ - t  05  sm4<£  -  30  am £tp  +  3) 

R  I  r  1  r3  \  3  /  35  r5  J 


where:  G  =  universal  gravltion  constant  =  6.664  x  10"®  dyne  cmVgm2!  M  =  mass  of  the  earth  « 
5.975  x  1027  gm;  R  a  equatorial  radius  of  the  earth  =  6.3783  x  10®  cm)  Kj  >»  1.637  X  10  ®;  and 
K2  »  1-07  x  10”®.  The  force  of  gravity  is  then  given  hy 


+■ '  ‘  -  [r»  F 7  &  ~ *  5] 


7 m  8 1 (r . <£)  ~  ?o «2 (r > & 


8l(r,<f> )-  91\JL+3 K,  *2.  ^i(35sin4^-30sin2^+ 3M  (20) 

R  Lf2  r4  v  7  f6  J 

g2(r,^)  ■  -  2K2  ^2  sin  0  cos  ^  (21) 


The  quantities  Xj  and  arc  a  measure  of  the  earth’s  oblateness.  When  Kj  =  A!2  =>  0, 
f =  -T0  (GAf/r2). 

Equation  (4)  can  be  written 


B+D+L+N+J+  mg  +  /q  =  m  ^  —  +  20x7+0x0x7"^ 


“Jeffreys,  H.,  The  Earth:  Its  Origin,  History  and  Physical  Constitution,  3rd  ed.,  p.  129.  University 
Press,  Cambridge,  England,  1952. 
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Substituting  Eqs,  (14),  (15),  (16),  (17),  and  (19)  in  Eq.  (22)  yields  the  following  three  Equations: 


--  -  B  -  ])  sin  (6  +  a  +  17)  -  D  sin  6  +  L  cos  ®  -  N  cos  (©  +  a)  -  mg  j  j 


v  cos  ®  0  +  v  sin  ®  -  —  cos2  ( 
r 


+  2fiv  cos  ®  cos  4>  sin  0  -  fi2r  co»^  <f> 

'  ■  \ 

..V, 

~[{/0-fl-7)  cos  (8  +  a  +  t?)  cos  a  -  D  cos  8  cos  a  -  L  sin  ®  cos  a 

+  IV  sin  (0  +  a)  cos  a  -  mg2j 
a  ~  v  cos  8  sin  a  b  -  v  cos  a  sin  8  8  +  v  cos  ®  cos  a 


1^2  ^2  9  o 

+  sis  8  cos  8  cos  o  +  —  cos2  8  sin2 17  tan  </> 

T  T 


-  2  fl  v  cos  8  sin  qi  Bin  a  +  02r  sin  t/>  cos  <f> 

—  |(/0  -  B  -  /)  cos  (8  +  a  +  tj)  sin  a  -  D  coo  8  sin  a  -  L  sin  8  sin  o 
n  L 

.  +  N  sin  (8  +  a)  sin  oj 

o  v  cos  ©  cos  a  a  -  v  sin  a  sin  8  8  +  v  cos  ©  sin  a 


2  2 

+  ~  sin  8  cos  8  sin  0  -  —  cos2  ©  cos  o  sin  0  tan  <j> 
r  r 


+  2  Civ  cos  8  sin  </>  cob  a  ~  2(1  v  sin  ©  cos  <j> 
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Solving  the  preceding  hree  Equations  for  v,  0,  p  gives 


”  a  —  £(/0  -  B  -  J)  cos  (a  +  ij)  -  oj  +  ~  sin  a  -  g  j  (r ,  <A)  sin  0  -  gjj  (r ,  ^i)  cos  o.  cos 
-  O^r  sin  2<£  cos  0  cos  a  -  sin  0  cos®<£^ 


© 


(23) 


© 


]  |'  T  /  t,2  \ 

—  |y0  -  B  -  /)  sin  (a  +  ij)  +  L  -  N  cos  aj - f 1  -  - — \ 


g2 


02-  /  1  . 

- co®  a  ®in  ®  -  2Q  co®  <j>  sin  a  +  (  —  sin  26  sin  ®  co®  cr 

v  v  \  2 


+  cos^  cos  0  j 


(24) 


a  *  —  cos  0  sin  oiontj>~  20  (sin  tf>  -  tan  0  cos  <j>  cos  a) 


~  £*2  +  y  O2''  sin  2^J  sec  0  sin  o 


(25) 


Let  M  be  the  sum  of  the  torques  tending  to  rotate  the  missile  about  a  transverse  axis 
passing  through  its  C.  M.,  normal  to  Tend  r,  and  let  /  be  the  moment  of  inertia  about  this  axis. 
Then,  from  Eq.  (5), 


--[/(©+  a)3  -  m  (0  +  a)  =  M 


(26) 


Further,  let  T  be  the  sum  of  the  torques  (due  to  a  ring  of  jets  placed  around  the  circumference  of 
the  missile)  tending  to  rotate  the  missile  about  its  longitudinal  axis,  and  let  /'  be  the  moment  of 
inertia  and  <u  the  angular  velocity  about  this  axis.  Then,  from  Eq.  (8), 
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—  (/ '  ai)  -  ma  ( ajo  -  v  )  =  T 

dt  a 


(27) 


where  a  is  the  perpendicular  distance  from  the  longitudinal  axis  to  a  jet  and  t>g  is  the  exhaust 
velocity  of  the  jets. 

Equations  (13),  (23),  (24),  (25),  (26),  and  (27),  plus  an  appropriate  control  (guidance) 
equation  in  pitch,  completely  describe  the  motion  of  the  missile  relative  to  the  earth. 

If  (rQ,  <)>0,  0q)  la  the  launching  point,  and  (rj,  <f>j,  Oj)  is  the  impact  point  of  the  missile, 
then  the  impact  range  X/  is  given  by 


■  ■  cos“*  [  cos  <f>Q  cos  coa  (0q  -  6j)  +  aln  tf>Q  sin  <j>j] 


(28) 


Here,  it  is  assumed  that  rj  ■  r0 . 

B.  Satellite 

The  missile,  or  e  part  of  it,  may  permanently  (Ideally)  move  about  the  earth  in  a  nearly*5 
circular  orbit,  provided  that  v  and  0  assume  appropriate  values  at  the  instant  of  final  thrust 
termination.  These  values  are  determined  as  follows: 

If  tjie  orbit  were  perfectly  circular,  then  r(t)  a  0,  for  t>  t\  where  t'»  time  of  final  thrust 
termination.  It  then  follows  from  Eq.  (18a)  that 


0(1)  b0  (!>!') 


..  /  - 

and,  consequently, 


(a) 


0 (t)  a  0  UHt*) 


(b) 


Applying  conditions  a  and  b  to  Eq.  (24)  yields 


v  =  fir  cos  <f>  sin  a  +  (gjr  -  fl2r2  cos2  <j>  cos2o)*(  (l  >  1 ') 


'’The  oblateness  of  the  earth  prevents  the  orbit  from  being  perfectly  circular  and  also  prevents  it 
bom  lying  in  a  plane. 


-i/ 
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which  is  the  speed  of  the  satellite  relative  to  the  earth.  Thus,  if  the  satellite  is  to  desbribe  a 
nearly  circular  orbit  at  a  distance  ■vr  from  the  earth's  center,  then 


v(t ')  -  Or(t ')  cos  ^ (< ')  sin  o  (< ')  +  [gx ( t  l)r(t ')  -  0 V  (< ')  cos ')  cos 2o  (t  ')]*  (29a) 


©U')  =  0 


(29  b) 


Equations  (19),  (23),  (24),  and  (25),  with  fa  -  fl  ■  /  ■  L  »  N  •  0,  therefore  deaorlbe  the  motion  of 
the  satellite  relative  to  the  earth,  subject  to  the  "initial"  conditions 


r(t') 

tl>W) 

eu1) 

Eq.  (29a) 
Eq.  (29b) 
o(t‘) 


V  (30) 


If  r(t')  is  large  (e.g,,  r(l')  -  1,25  R),  then 


«1  -  2 


82  -  0 


D  -  0 
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In  such  a  case,  the  orbit  could  be  considered  circular,  so  that  the  equations  of  motion  would  be 


v  sin  o 


r(l')  cos  4> 


a  m  ^JL-.  sin  a  tan  d>  -  2fl  sin  0  +  ^  ein  20  nib  0 


rU‘) 


z  v 


>.  (ai) 


where 


v  -  Or(t ')  non  0  sin  <t  + 


-  0s  ra(t')  con  2  0  ooa2ff| 


(32) 


The  value*  of  0(t')  and  a(t')  determine  the  initial  angle  that  the  orbit  mokea  with  the 
equatorial  plane.  Thus,  in  Fig.  8,  the  orbital  angular-velocity  vector  "to  of  the  satellite  makes  on 
angle  0  with  the  rotational  aula  of  the  earthi  X  is  the  angle  betweenvthe  positive  x  axie  and  the 
projection  of  a  in  the  x-t  plane. 


Then 


^  am 

a  «  ka  sin  0  cos  X  +  jeo  cos  0  +  ka  ein  0  sin  X 


Initially,  T  lies  in  the  plane  of  the  orbit,  so  that 


0>  =  1"3»|  =  -y  = 


In  accordance  with  Eq.  (11), 


a  =  70  6>r  +  0O  W(^  +  T0 
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where 


o>r  =  a)  bid  ip  coa  <p  coa  (A  -  6)  +  oj  cos  <p  sin 


(33a) 


a>^  »  -  fi)  sin  ip  sin  <p  cos  (A  -  0)  +  <y  cos  ip  coa  <p 


(33h) 


<u  0  =  m  ein  ip  ain  (A  -  0) 


Now, 


(f-Sx? 

it 


(33c) 


where  d/rft  refers  to  inertial  apace.  If  the  earth  is  takes  as  the  reference,  the  equation  above 
becomes 


dT  ~  ^  _  * 

—  +  u  x  r  *■  <u  x  r 


But,  from  Eq.  (12), 


70t+<P0v  cos  0  cos  a  +  v  cos  0  sin  o 


Further,  0,  x’T=  -  0oOr  cos  <?!>,  and  «  x7-  4>Qr<o6- e0r<Of  Substituting  in  Eq.  (34)  and  equating 
components  yields 


T  =  0 


—  cos  0  cos  a 
r 


=  (1  cos  <p  -  —  cos  0  sin  o 
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However,  for  t  ~  I 8  «  0.  Combining  the  last  two  Equates  with  Eqa.  (83b)  and  (33c)  gives 


—  cob  a  =  (i)  sin  sin  (X  —  0) 
r 


(35a) 


0  cos  -  >—  sin  o  =  o)  cos  \ji  cos  <f>  -  0)  sin  ^  sin  <f>  cos  (X  -  6) 


(35b) 


Sines  Jo  •  T®  0  initially,  o>r  »  0.  Then,  from  Eq.  (83a), 


cos  (A  sTh  “  -  sin  co*  <f>  cob  (X  -  d) 


(35c) 


cos  ^  j  ^  sin  o  cos  qi  +  (1  cos2q5j 


Consequently, 


cos  (fr(t')  «*  H{t‘) 

--N(t') 


(-ir<iff(t')  <  0) 

(0  <  o  (t ')  <  ff) 


whore 


H(t')  - 

*!<•) 


.7.5^.—  sin  ff(t ')  cos  ^(t*)  +  0  eos2^(t') 
r(t') 


’] 


Here,  v(t')  is  given  by  Eq.  (29a). 
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VI.  APPROXIMATE  LOCATION  OF  THE  MISSILE  IMPACT  POINT 
FOR  A  ROTATING  EARTH 

The  exact  location  of  the  miaailo  impact  point  ia  obtained  by  integrating  Eqa.  (13),  (23) 
to  (27),  and  the  appropriate,  control  equation  from  takeoff  ( t  =  t^)  to  imnaot  (t  =  </).  It  ia  poasible 
to  deacribe  the  approximate  location  of  this  point  by  integrating  the  equations  only  up  to  t  =  t  , 
the  time  at  which  the  guidance  ia  terminated  (the  thrust  hewing  terminated  at  t  5  *')■  The  impact 
point  can  then  be  determined  from  the  values  of  the  missile  position  and  velocity  at  t  ■  t  i 
«•(<'),  v(t'),  ©(«')#  and  ff(t'). 

It  ia  assumed  that  the  earth  ia  a  perfect  sphere  (inverse*sqnare  force  law),  and  that  the 

mlaSilc  mu V c a  through  a  Vauuuut  during  the  iuiurVai  t  5  ft  —  £j.  llctiCc,  it  W ill  diownuc  Su  Cirippu 

(in  inertial  space)  during  this  time  (aee  Fig.  9).  Here,  /j  and  /2  are  the  foci  of  the  ellipse,  which 
lies  in  .a  plane  fixed  in  Inertial  space.  Let  $  be  the  angular  displacement  of  the  missile  from  the 
line  passing  through /j  and /j.  Then  the  equation  of  the  ellipse  is 


JL  -  1 

r  s 


( 1  -  e  cos  $) 


(36) 


where 


o^VVJcos2©^') 


am 


80r0 


92*  1 


2* 


r(l')  r2(t')  cos20i(t') 


The  subscript  t  denotes  that  the  associated  quantity  is  measured  in  inertial  space.  Here,  e  is  the 
eccentricity  of  the  ellipse,  and  g0  and  are  the  values  of  g  and  r  measured  at  the  earth's  surface. 
(The  impact  altitude  is  assumed  equal  to  die  launching  altitude.) 

The  impact  point  <f>(Uj),  is  a  function  of  r(t'),  fyft'),  u((t'),  ©^U'), 

ff^(t')  (refer  to  the  spherical  right  triangle  in  Fig.  10).  From  Fig.  10, 


cos  o({t ')  = 

tan  x 


(37a) 
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■in  Ojit ')  m  — -  —  (37b) 

«n  X 

Substituting  this  in  Eq,  (37b)  yields 


cos  £  m  cos  ^(W)  oos  $y(ly)  cos  [  9^t')  -  fy(t/)]  +  sin  ^(t ')  sin 


From  Eq.  (28), 


oos  oos  tftfUf)  ooaffyO')  -  OtUj))  +  tin  ^00  sin  flbjU;)  -{1-  sinaff<(t')sio2x]',< 


From  Eq.  (28), 


Combining  Eqs,  (38)  and  (89)  gives 


(88) 


cos  oos  $)/(*/)  cos  [ 9 fit *)  -  flj(i/)l  +  sin  qSy(r')  sin  <fy(t/)  -  cosy  (89) 


sin  #{(t() 


( 1  -  sin2ffyO ')  sin2*)  cos  ^j(t ')  -  cos  $t(t ')  cos  * 
sin  C  *,(*')  -  *,(*')! 


(40) 


But  17  =  ^y(t')  -  ^(<0.  Hence,  from  Eq.  (87s), 


}^y(t ')  b  ^y(t')  +  ten'1  [  cootf/  (<')  ten  *] 


Substituting  this  in  Eq.  (40)  yields 


sin  <£y(t/)  =  cos  ^(t')  cos  Oy(t')  sin  *  +  *in  4>jW)  cos  * 


(41) 
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ti|  =  ii  t  Qxr 


•T?"  ) 


Now,  fl  =*q  fl  sin  <f>  +  ^  008  ^  •  Substituting  this  expression  and  Eq.  (10)  into  Eq.  (44)  yields 

the  three  Equations 


v(  sin  0^  =  t>  sin  0 
v(  cos  cos  Of  «>  v  cos  6  oos  o 


Vj  COB  Oj  BIO  Of  m  V  OOS  O  sin  u-  —  f  y- 


which  combine  to  give 


■’w  JL 


vj  -  vi  +  Or  cos  0  (Or  cos  qi  -  2v  cos  0  sin  a) 


From  Fig.  11, 


sin  0j  ■  sin  0 

*  li  _ 


tan  a i  «  tan  o  -  0  —  sec  0  sec  a  cos  <jt 


x(  -  *  cos  0  (<-  l1)  +  s  sin  ft  (t  -  t1) 


1 1  •>-  x  sin  fl  (t  - 1 ')  +  *  cos  fl  (t  ■■  t ') 


Substituting  Eqs.  (9)  into  the  <  :pressions  above  yields  the  three  Equations 


cos  9 1  cos  q&j  =  cos  9  cos  <j>  cos  fl  (t  -  t ')  +  sin  9  cos  <j>  sin  fl  (t  -  t ') 
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sin  tf>i  *»  sin  ^ 


sin  6,  cos  d>i-~  bob  6  cos  0  sin  0  (*  -  *  0  +  »iu  <5  cos  </>  cos  fi  (t  -  t') 


which  combine  to  give 


dt-  e-n(t-t') 


Equations  (45)  and  (46)  ore  the  desired  relations.  The  quantity  t  -  t'  is  obtained  as  followsi 
from  Eq.  (36),  . 


-V-  —  sin  <&  •  4 

r  * 


Combining  this  and  Eq.  (86)  with  the  Equations 


v{  sin  -  r 


v{  oos  @|  -  r  4> 


p  m  rpj  com  ©^  *  angular  momentum  ■  constant 


h.  --+  ^[(e2-!)  4  +— "  l]  U*S.Ktt 

—  it  r  I  a2  4  j 


—  f(«2  -  1)  — +  —  -  1 

r  L  s2  s 


(».,  <  t  <  */)  (47b) 
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Here,  t$  =  time  at  which  the  missile  reaches  the  summit  (see  Fig.  9).  Integrating  Eq.  (47a)  from 
t  to  <4  and  Eq.  (47b)  from  t4  to  t  (t4  <  C  £  tj)  and  combining  the  two  results  gives 


t  -  t' 


H 

4  Jr(t) 

i  1  h  M 

+  r(t') 

2  _ 

l  4  \ 

ib-){  1 

e2  V  r(t)  )  _ 

~  r(t')) 

+  — i_.«n -l[jL(e-±)r(t)+JL] 

V— i—  (..JLU  *)  +  JL1 

aU-e8)*  i*  \  •! 


it*  | 

.(!-.¥  J 


(48) 


.Hera,  Pm  rU  ')vt(t')  eoo  6j(t . 

Equations  (43),  (48),  (46),  and  (48),  used  la  conjunction  with  Eqa.  (41)  and  (42),  yield  the 
desired  impact  point  ^(t;),  0(t/). 

H,  for  simplicity,  it  la  assumed  that  the  earth  does  not  rotate  daring  the  Interval 
*0  -  *  -  *  •  then  the  miselle  path,  during  this  Interval,  will  He  in  a  plane  relative  to  the  earth) 
therefore,  the  equations  to  be  integrated  will  simply  bei  Eq.  (18a)|  the  expresaion 
n{Wr)v  008  ®  (where  *  -  ground  range  measured  along  the  adrfaee  of  the  earth))  Eqe.  (23)  and 
(24)  (with  0  ■  0))  Eq.  (26))  and  a  control  equation  in  pitch..  However,  only  the  quantities  r(t (), 
v(t')  and  8(t')  are  obtainable  from  the  solutions.  The  remaining  quantities  <t>W),  d(t'),  und 
er(t')  can  be  determined  from  the  known  values  of  ^(*0),  0(to),  o(tQ)f  and  *(t').  Thus,  analog 
to  Eqs.  (41)  and  (42), 


:ous 


<£(«')  ■=  sin-1  cos  0(to)  cos  cr(t0)  sin  sin  ^(t0)  coo 


’r0 


r0 


(49) 
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&(t  )  =»  9(t0)  +  co*"1  |sec  ^(t0)  sec  ^&(t‘) 


*(t0 


-  tan  4>  (<0)  tan  #  (<  *) 


'3 


^(*q)  -  W*“  [  ] 

*<‘o>  ! 


(0  <  a  (t0)  <  77) 


(“7 T<o(t0)  <  0) 


(o  (t0)  m  0,  77,  -  77)  (50) 


V 


■\ 


Consider  the  spherical  right  triangles  In  Fig.  12.  From  this  Figure, 

tan  (3j 


tan  o(tg) 


sln^(t0) 


tan  ff(t') 


tan  (3  2 
ela$(i ') 


tan^(*0)  tan^(f') 


aln/3j  alnjBj 


These  combine  to  give 


0 (*  )  *  *ln  1 C  cos  ^ (i0)  tin  0  (*0)  aeo  </>(t ')! 


(51) 


Th.  .nll«  ft.  more  .Ill  ft.  tap..,  p.tatftf,,),  «(,,).  ,.  .kftto.d  (torn  ft... 

..1...  of  r(, '),  <ft(t «(,'),  «(.-).  8(,'),  ud  ,.f..M.  „kk  ft.*,,,), 

In  the  preceding  paragraph. 
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Table  1.  Nomenclature 


a 

C.  M. 


F  - 
flm 
to  " 
f\t  f 3  " 

G  ■ 

7  - 

— *0-  - 

/  - 

<»  /.T  - 
7 . 

n. 

*  J 

mdm 

L  - 
t  a 

M  - 

7  = 


perpendicular  distance  from  the  missile  longitudinal  a  as  to  a  jet  located  on  the 
circumference. 

base  drag  force. 

center  of  mass  of  the  missile. 

velocity  of  mj  with  respect  to  the  C.  M.  of  the  missile, 
drag  force. 

eccentricity  of  ellipse  (Fig.  9). 

external  forces  acting  on  mj  and  m. 

thrust  from  the  fth  jet. 

vacuum  thrust. 

fool  of  ellipse  (Fig.  9). 

universal  gravitation  constant. 

acceleration  of  gravity. 

value  of  gravity  acceleration  at  the  launching  point.  _  . 

moment  of  Inertia  of  the  missile. 

unit  vectors  pointing  in  the  direction  of  increasing  x,  y,  *,  respectively, 
jet  vane  and/or  exhaust-obstruction  effect. 

'll 

constants  which  are  a  measure  of  the  earth's  obleteness. 

angular-momentum  vector  (Sec,  III), 
lift  force  (Sec.  V). 
mass  of  the  earth. 

torques  taken  about  a  transverse  axis  through  the  C.  M.  of  the  missile, 
mass  of  the  missile. 
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Table  1  (Cont'd) 

«  mass  of  a  quantity  of  exhaust  expelled  from  the  missile. 

=  total  number  of  jets, 
a  lift  force  of  the  rudder, 
a  position. 

=  orbital  angular  momentum  -  constant. 

«  total  linear  momentum  of  mj  and  m,  measured  in  £  (x,  y,  *). 
a  equatorial  radius  of  the  earth. 


a  position  veotor  drawn  from  the  center  of  the  earth  to  the  C.  M.  of  the  missile. 

a  impact  point  of  the  miasile. 

>  launching  point  of  the  missile. 

a  unit  vectors. 

a  unit  vectors  pointing  in  th  '  ection  of  increasing  r,  <j>,  0,  respectively, 
a  a  constant  determined  by  the  missile  position  and  velocity  at  t  =  t ' . 

“  average  distance  from  the  C.  Vi.  of  the  missile  to  the  base  of  the  missile. 

=  position  vectors. 
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labie  1  (Conf'd) 


•i  >  =  unit  vectors 


t  =  time. 


t '  m  time  of  final  thrust  termination  and/or  time  of  pitch-guidance  termination. 
tf  x  impact  time  of  thn  missile. 

*  time  at  which  the  missile  reaches  the  summit  of  its  irajectory. 
t0  »  instant  of  time  Immediately  before  m j  la  expelled  (Secs  II). 
t0  +  A  »  Instant  of  time  immediately  after  mj  is  expelled  (Sec.  II). 

(q  m  launching  time  of  the  missile  (Sec.  VI). 

V  »  gravity  potential  function  for  an  oblate-spheroid  earth. 

v  -  |T|. 

T (()  x  velocity  of  the  missile  C.  M.  -  dr7di. 

ve  »  exhaust  velocity  of  jets  situated  around  the  circumference  of  the  missile. 
X  x  ground  range  measured  along  the  surface  of  the  earth. 
x /  x  impact  ground  range  of  the  missile. 
xi>  y i'  zi  “  Cartesian  coordinates  measured  in  inertial  space, 
a  x  angle  of  attack. 

f*i[ 

>  *  arcs  of  great  circles  on  the  earth's  surface  (Fig.  12). 

P  2J 

y  =  angular-velocity  vector. 

jj  =  angle  between  the  motor  thrust  and  the  missile  axis  (Sec.  V). 


i 

i 
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Table  1  (Cant’d) 


»? 

* 

C 

*< 

8  « 
0  - 

X  - 

X  U,  y,  *)  *> 
S '  U  y ', « ')  - 

a  “ 

T  ■= 

0  - 

<t>r ] 

«i 

vt 


<t>  C9 

t/y  t* 

a  - 


O)  = 

6)  ss 


arcs  of  great  circles  on  the  earth's  surface  (Sec.  VI). 


inclination  of  7  to  the  local  horizontal, 
longitude. 

angle  between  the  +  x  axis  and  the  projection  of  7  in  the  xz  plane  (Fig.  8), 

set  of  orthogonal  axes  fixed  In  inertial  space. 

set  of  orthogonal  axes  fixed  in  the  earth. 

angle  between  the  local  horizontal  and  the  local  meridian. 

torques  taken  about  the  longitudinal  axis  of  the  missile. 

latitude. 

quantities  measured  in  Inertial  space. 

angle  between  rand  the  line  through  /j  and  /j, 

angle  between  a  and  Cl  (Fig,  8). 

spin  angular-velocity  vector  of  the  earth. 

angular-velocity  vector  of  the  miBBiie  about  the  longitudinal  axis  (Sec.  III), 
orbital  angular  velocity  of  the  satellite  (Sec.  V-B  and  Pig.  8). 
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Fig.  10.  •  Missile  Impact  Point  «i  j. 
Projection  of  Missile  Shutoff 
Point  on  Earth’s  Surface 


CENTER  OF 


Fig.  11.  Rotation  of  the  Earth-Fixed  System 
,  Relative  to  Inertial  Space 


V  (/') 

v 


•  EQUATOR 


Fig,  12,  Spherical  Triangles  on  Earth’s  Surface  Relating 
o  (i  ’)  and  tr  (t„) 


